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ABSTRACT 

Self-gravitating accretion disks collapse to star-forming(SF) regions extending to the inner edge of the dusty 
torus in active galactic nuclei (AGNs). A full set of equations including feedback of star formation is given 
to describe the dynamics of the regions. We explore the role of supernovae explosion (SNexp), acting to 
excite turbulent viscosity, in the transportation of angular momentum in the regions within lpc scale. We 
find that accretion disks with typical rates in AGNs can be driven by SNexp in the regions and metals are 
produced spontaneously. The present model predicts a metallicity-luminosity relationship consistent with that 
observed in AGNs. As relics of SF regions, a ring (or belt) consisting of old stars remains for every episode 
of supermassive black hole activity. We suggest that multiple stellar rings with random directions interact and 
form a nuclear star cluster after episodes driven by star formation. 
Subject headings: accretion disks - quasars: general - stars: formation 



1. INTRODUCTION 

Accretion onto supermassive black holes (SMBHs) is the 
energy house of active galactic nuclei (AGNs). Applying 
models of accretion disks to fit the observed big blue bumps 
in the continuum of AGNs and quasars, one finds that SMBHs 
accrete typically with sub-Eddington rates (Laor & Netzer 
1989; Sun & Malkan 1989; Brunner et al. 1997; Lii 2008). 
Observational examination of the Soltan argument for the cos- 
mological growth of SMBHs prefers that they grow episodi- 
cally through baryon accretion (Yu & Tremaine 2002; Wang 
et al. 2009a), in agreement with the paradigm of coevolution 
with their host galaxies. However, what is governing accre- 
tion rates of active SMBHs? 

We have several relevant clues to understanding the for- 
mation of the accretion flows approaching SMBHs at ~ lpc 
scale. First, the well-known metallicity and luminosity rela- 
tionship in AGNs and quasars (see a review of Hamann & Fer- 
land 1999) strongly implies that the flows may undergo fast 
metal productions through stellar evolution in active episodes 
in light of the absence of cosmic evolution of metallicity (Di- 
etrich et al. 2003; Shemmer et al. 2004). This gets sup- 
ports from evidence that there is intensive star formation with 
the top-heavy initial mass function (IMF) within ~ O.lpc in 
the Galactic center (Paumard et al. 2006; Nayakshin & Sun- 
yaev 2005) confirmed by numerical simulations (Nayakshin et 
al. 2007; Bonnell & Rice 2008; Hobbs & Nayakshin 2009). 
We may naturally postulate that accretion onto SMBHs is the 
consequence of metal production. Second, it has long been 
known that accretion disks in quasars are massive and self- 
gravitating in outer regions (Paczynski 1978; Kolykhalov & 
Sunyaev 1980; Shlosman & Begelman 1987; Collin & Zahn 
1999). Obviously, the self-gravity is driving intensive star for- 
mation ongoing there (Goodman 2003), and a star-forming 
(SF) disk (or regions) appears. It seems that there is no way to 
avoid the production of metals and trigger accretion flows in 
the SF disk, yielding the metallicity-luminosity relationship. 
We argue that SNexp plays a leading role in the establishment 
of the relationship in light of SNexp-excited turbulent viscos- 
ity as shown in the ~ lkpc regions by numerical simulations 
(Wada & Norman 2002; Hobbs et al. 2010) or physical moti- 



vation (Wang et al. 2009b; see also Kumar & Johnson 2010) 
and - 10 4 Sloan type 2 AGNs (Chen et al. 2009). Further- 
more, we are inspired by the recent evidence for that the "vis- 
ible" inflows displayed by the H/3 line are likely developed 
from the inner edge of the dusty torus in ~ 5000 Sloan quasars 
(Hu et al. 2008a,2008b). All clues merge into a question: are 
accretion inflows driven by star formation in the massive self- 
gravitating disk? 

In this Letter, we show the feasibility that an SF disk is able 
to drive a Shakura-Sunyaev disk (hereafter SS disk) in AGNs 
through SNexp-excited turbulent viscosity. We find that such 
a scenario agrees with observations of metallicity and results 
in the formation of a nuclear star cluster potentially. 
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FIG. 1. — Schematic illustration of the present model. The SF disk starts 
from the inner edge of the dusty torus (R sl ,b) and switches on the SS disk at 
the self-gravity radius (Rsg). 

2. BASIC EQUATIONS 

Figure 1 shows an illustration of the present model. The 
outer boundary of the SF disk (R out ) is chosen at the inner 
edge of the torus as suggested by Hu et al. (2008a,2008b), 



where dust particles are sublimated at R su ^ = 1.3L 
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where Luv.46 =iuv/10 46 erg s" 1 (Barvainis 1987). The op- 
tion of R out = /? su b is motivated by a fact that emission from 
R > ^sub is invisible in optics. Collisions of molecular clouds 
transfer angular momentum outward and supply gas to the 
SF disk with a rate of M out = 1 ~ 20M Q yr _1 at R su \, from the 



torus (Krolik & Begelman 1988) with a typical mass range of 
10 5 — 1O 7 M (Mor et al. 2009). This typical mass is just at 
the level required by a single episode of SMBH activity. 

2.1. SNexp-driven gaseous disk 

With the inclusion of the mass dropout from star formation 
and the injection of gas recycled from SNexp, the mass con- 
servation equation for the SF disk read£] 
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where S gas is the surface density, Q = — 27ri? 3 i/£g as (dVt/dR) 
is the viscosity torque, v is the kinematic viscosity, ft is 
the angular velocity, E* is the surface density of star for- 
mation rates, and Esn is the injection rate of gas recycled 
from SNexp. We use the relation Esn = /r^X, where f c is 
the fraction of the SNexp-ejected mass to the total formed 
stars. For an episodic activity of AGN with A?g (e.g., Wang 
et al. 2006), the minimum mass of the stars producing an 
SNexp should be M C /M Q > 7.0 Af^] 4 , where the lifetime of 

hydrogen main-sequence stars is ?ms = 13 (M*/M Q ) ' Gyr 
and Af .i = Af G /0.1Gyr. Assuming an IMF as N(M*)dM* = 
NQhT/dM^, we have f c = J^ M@ M^N{M^)dM^/M tol , where 

Mot = J M ° M*N(M*)dM* is the total mass of the formed 
stars. It follows f c ~ 0.38 for the Salpeter IMF (3 = 2.35, 
M m i n = 1.0M Q , and M c = 7M Q , meaning that a significant 
mass fraction of the formed stars will be injected into the 
SF disk from SNexp (here we neglect the remnant compact 
stars after SNexp). The scale influenced by an SNexp is of 

the order of Dsn ~ {Esn/ni p ncf) ~ 0.1 E 5l (nnT^Y - ' 
pc, where £51 = £sN/10 5I erg is the kinetic energy of SNexp, 
«n = Hgas/lC'cm" 3 is the number density of the medium in 

the SF disk, c, w 3.0 x 10 5 r 3 cm s" 1 is the isothermal sound 
speed, and T^ = T/10 3 K is the temperature. We find from 
this simple estimation that the length of the SNexp influence 
is comparable with the height of the SF disk (~ 0.1 pc) and 
the SNexp-excited turbulence velocity V tur 3> c s , and hence 
SNexp plays a key role in the transportation of angular mo- 
mentum. 

In this Letter, we focus on the stationary case, namely, 
<9E gas /<9f = 0. Introducing the parameter X = R 2 il, Equation 
(1) is rewritten as 
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(2) 



where the Kennicutt-Schmidt (KS) law S* = ASg as is used for 
simplicity in the entire SF disk, A is a constant, and 7 is the 
index (Kennicutt 19980. Lynden-Bell & Pringle (1974) show 

Equation (1) is the modified version from Lin & Pringle (1987) and Wang 
et al. (2009b). We point out that Ssn is at f— T*, where r* is the hydrogen 
main-sequence lifetime of stars depending on the stellar mass. Since stars 
form in the same gaseous disk, the stellar disk always follows the latter in the 
current model. Unless the stellar disk was formed in the previous episode, 
it could be larger than the gaseous disk, leading to reduce the feedback effi- 
ciency significantly (Nayakshin et al. 2007). 

2 We note that star formation rates may deviate from the KS law along the 
radii in the SF disk. It is not sufficiently understood (Krumholz et al. 2009), 
especially for so dense disks. 



that dQ/dX = M at any radius, where M is the mass rate of the 
inflows, provides the inner and outer boundary conditions. We 
employ the a-prescription viscosity as ^ SN = a SN V tnr H, where 
a SN is a constant, V tul - is the turbulence velocity driven by 
SNexp, and H is the thickness of the disk (Wada & Norman 
2002; Wang et al. 2009b). The energy equation of the SNexp 
excited turbulence is given by 
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— — — = e/*p*£sN, 



(3) 



where fdis = H/V tur is the timescale of the turbulence, e is the 
efficiency of converting kinetic energy of SNexp into turbu- 
lence, and /*/9* is the density of SNexp rates. Here /* is a 
number fraction of stars to the total mass of the formed stars, 
which are able to produce SNexp during one AGN episode. 

For a given IMF, we have /* = J M ° N(M*)dM*/Mtat and 

/* - 2.3 x IO^Mq 1 for j3 = 2.35. Considering S* = p\H and 
Sgas = Pga&H , we have 



^ = e /,A£s N S^ 1 - 



(4) 



Following Paczynski (1978), we assume that gas rotates 

1 il 

with a Keplerian velocity J7k = (GAf.//? 3 ) , where G is 
the gravity constant and M m is the SMBH mass. A ver- 
tical equilibrium holds among the disk self-gravity, SMBH 
gravity, and turbulent pressure. We usually have dP lm /dz = 
(fi 2 z+27rG£ga S )/9gas, where P tm = p g a S V t u r is the turbulence 
pressure and the second term is the self-gravity. Since the 
star formation supports the thickness of the SF disk, the self- 
gravity can be neglected (Thompson et al. 2005), allowing us 
to have an approximation form of vertical -averaged structure, 
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We point out that the above equations are actually averaged 
in the vertical direction, and the parameters obtained in these 
equations are the values at the mid-plane of the SF disk. 

The Q - S g a. s relation is specified for further simplification 
of Equation (2) through the viscosity v. With the help of 
Equations (3)-(5), we have 



dX 2 



B 

X 2 



= 0, 



yielding an analytical solution in the form of 

Q(X) = Cl Xexp f I J +c 2 Xexp f~ 



(6) 



(7) 



where B = 2GM. [(l-f c )/3a SN ef*E SN ] 1/2 , X = (GM.R) 1/2 
for the Keplerian rotation, and c\ and C2 are two constants de- 
termined by the boundary conditions. At the outer radius of 
the SF disk, only the mass injection rate is known and given 
by the mechanism of molecular cloud collisions (Krolik & 
Begelman 1988). We have dQ/dX\ x =M out , namely, 



d (1 -BX£) e^ +c 2 (1 +BX^) e~* 



■■M n 



(8) 



1 11 
where X oat = (GM m R out ) ' . We assume that the inner radius 

of the SF disk is set to be the self-gravity radius (R sg ) of the 
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SS disk, which depends on the SMBH accretion rates (M.). 
This yields dQ/dX\ x =M,, and subsequently 

c\ (l -5X" 1 ) e^i + c 2 (1 +BX^) e' 1 ^ = M„ (9) 

where X in = X sg = (GM,R sg ) I/2 . 



rates of the SS disk as well as the structure of the SF disk are 
obtained self-consistently. 




FIG. 2. — Structure of SF disks, which drive SS disks with an accretion 
rate of 0.5 — 2Lrj<id/ c . as shown by the solid lines in ay. The panels on the 
left column are for e = 0.5, but M out takes 1.0, 2.0, 3.0, 4.0, 5.OM0vr' repre- 
sented by the black, red, green, blue, and cyan lines, respectively. The dashed 
lines are the local star formation rates defined by M, = 7ri? 2 S*, which are 
of 1 - lOMQyr -1 . The panels on the right column are the structures for fixed 
Mont = 3.0M o yr', but e = 0.2, 0.3, 0.4, 0.5, 0.6 by black, red, green, blue, 
and cyan lines, respectively. 

2.2. Switching on the Shakura-Sunyaev disk 

The self-gravity instability develops when the Toomre's pa- 
rameter Q < 1, For the SS disk, the instability happens in the 
middle region of the SS disk at the radius 
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where /?sch = 2GM,/c 2 is the Schwartzchild radius, c 
is the light speed, M% = M./10 8 M Q , m = M m c 2 /L E dd = 
M,/0.2Mft,MQyr~ l , Lead is the Eddington luminosity, and 
Oio.\ = ass/0-1 is the viscosity parameter in the SS disk (see 
Kato et al. 1998). We assume that the SF process quenches 
within R sg , and the SF disks smoothly switch on the SS disks 
at R sg in term of surface density. The surface density of the 
SS disks is 



Ess = 8.2 x 10 8 cT 4/5 Mg /5 m 3/5 r 3 3/5 M Q pc" 



(11) 



where rj, = R/10 3 Rs c i 1 , giving rise to the expression of the con- 
dition (Sgas = Ess) of switching from the SF regions to the SS 
disk 

Cl e^i+c 2 e~^ = 37rae/*A£ S NSss^" 2 - (12) 

It is trivial to get the prolix expression of c\ and ci from Equa- 
tions (9) and (12), but we omit them here. We insert c\ and c<i 
into Equation (8) for the structure of the SF disk. Given R ont 
and M out , we can justify if the SF disk is able to drive the SS 
disk through solving equations (8), (9), (12). The accretion 
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FIG. 3. — Relation between black hole accretion rates and the star formation 
rates in the SF disks. Each line represents a relation for the labeled e. We find 
that a relation of 1Z t <x Af. and t determines intercepts of the relation. 

3. ACCRETION FLOWS APPROACHING TO SMBHS FROM SF DISKS 

The empirical KS law E* = AEg as is used, where 7 = 1.4, 
A = 2.5 x 10~ 4 , E* and E gas are in units of M yr -1 kpc~ and 
M pc" 2 , respectively. We employ £sn = 10 51 erg, a top-heavy 
IMF index = 1.9, f c = 0.63, /* = O.O3M \ ass = 0.2, and 
q;sn = 1 (in light of £>sn ~ H) throughout the paper. We take 
the outer radius R out = 1 .0 pc and M. = 10 8 M Q , but adjust two 
parameters e and M out to demonstrate the final accretion rates 
ofSMBHs. 

3.1. Structure of the SF disk 

Figure 2a 1 shows that the inflow driven by SNexp holds at a 
level of M. = 0.5 ~ 2M crit at the self-gravity radius (R sg ~ 0.02 
pc) for e = 0.5. This clearly demonstrates that an SF disk is 
able to drive an SS disk. Figure 2b\ plots E* along the radial 
direction. We find that there is a peak of E* in the SF disks 
before switching on an SS disk. This is caused by the star for- 
mation and the radial transportation of gas that compete along 
radii. From the second term of Equation (7), there is a cutoff 
of S + and then an SF disk switches on an SS disk. Figure 
2ci displays a smooth transition of the surface density from 
SF disks to SS disks. The bumps in Figure 2c\ are artificially 
made by the assumption for simplicity that there is no star for- 
mation in the SS disk. This assumption roughly holds since 
the temperature within R sg could be high enough to efficiently 
suppress star formation. 

The right column of Figure 2 shows a dependence of the 
structure on e for fixed M out . We find that M. increases with 
e, which appears in the exponential as shown in Equation (7). 
For M out = 3M Q yr _1 , the accretion rates are about 5%-10% 
M out from e = 0.2 - 0.6. The fraction ( 1 - e) of the SNexp ki- 
netic energy will be used to heat the gas of the SF disk, and 
then star formation is quenched close to the middle regions 
of SS disks. The output of this energy appears in the infrared 
bands in spectral energy distributions. Figure 3 shows the de- 
pendence of accretion rates on the global star formation rates, 
defined by 7?.* = J 2iril*RdR. For a fixed e, we find 7^* oc Mi, 
where q = 1 .43 - 1 .56 insensitive to M out . 
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FIG. 4. — Comparison of the present model with the observed metallicity- 
luminosity relation. We use the theoretical relation logZ/Z© = 1.0 + 
1.331og(N V/C IV) (Hamann et al. 2002) to convert N V/C IV into Z for the 
Shemmer et al. (2004) sample. We employ the relation Z4450 ~ 2L5100 ~ 
(2/9)Lbol, Aral ~ 9£<5I00 = riMtC 1 where r\ is the radiative efficiency and 
i-Bol, i-1450, and L5100 are bolometric and specific luminosities at 1450A and 
5100A, respectively. We use Z = 0, ey = 20, and e 2 = 0.1,0.05,0.01 for 
red, green, and blue solid lines, respectively, whereas the dotted lines are for 
£1 = 10. The theoretical metallicity-luminosity relation is consistent with the 
observed. 

It is necessary for AGNs and quasars to switch the SF disks 
from the middle region of the SS disks. We find that for lower 
M out the SF disks only support the outer region of the SS disk. 
In such a case, the SF disk may drive an advection-dominated 
accretion flow (ADAF) to SMBHs. Excellent examples of 
this case could be the Galactic center, some elliptical galaxies 
(e.g., M87), and LINERS, where star formation exhausts most 
of gas forming an ADAF as discussed by Tan & Blackman 
(2005). We would like to point out that (3 used here is only 
an intermediate top-heavy IMF. It has been suggested that (3 = 
0.85 (Maness et al. 2007) and /3 = 0.45 (Bartko et al. 2010) in 
the Galactic center. 

3.2. Relation Between metallicity and luminosity 

Metals are ejected through SNexp enhancing the metallic- 
ity in the SF disk. The mass of the metal element —i produced 

by SNexp is given by M Zi = f^dM*N(M*)ft*mf(M*,Z*), 

where m^(M^ : Z^,) is the ejected mass of the element—; de- 
pendent on the initial mass and metallicity of progenitor stars 
(Woosley & Weaver 1995; Gavilan et al. et al. 2005). 
With the Mz,, the total ejected metal is then calculated by 
Mz = z]2iMzr For simplicity, we take an approximation of the 
metal fraction m z = M z /M* w 0.1 -0.2 (Woosley & Weaver 
1986). 

Detailed calculations of metallicity involve evolution of 
stellar populations in the SF disk, advection of metal-enriched 
gas, and star formation, which make a strong radial depen- 
dence of metallicity complicated. However, we study the 
mean metallicity Z, as the observed, in the SF disk by in- 
troducing E\ and £2. Considering the radial dependence of 
metallicity, e\Z will be swallowed by SMBHs while ezZ is 
consumed by star formation. The lower limit of metallic- 
ity at radius R is roughly given by Z « f c m z A4*/M, and we 
have £i > Z in /Z ~ MfM/M*M t ~ M/M. ~ 5 from Fig- 
ure 2(ai), where the (sub/super) script "in" means the val- 



ues of parameters at the inner edge of the SF disk, and A4* 
is the local SF rates. On the other hand, its upper limit is 
Z max /Z ~ f c m z TZ*/M,Z ~ 30/(Z/Z Q ), where we use f c = 0.6 
and m z = 0.1. So we have 5 < e\ < 30 if the mean Z ~ Z Q . 
The parameter £2 ~ Zo/Z ~ 0. 1 for Z) = 0. 1Z . 

The total mass of the metal by time t through inte- 
grating the net increase of metals is given by Mz(t) = 

f (ZoM oal + f c 1Z*m z - £iZM, -etZIZ* J dr, where Zq is the 

initial metallicity of the gas supplied from the torus. The first 
term in the integration is contributed by the initial metal of 
supplied gas, the second is contributed by SNexp, the third is 
advected into SMBH and the last one is the consumption of 
metals during star formation. The total gas of the SF disk is 

given by M gas (f ) = f*" 2TrT, sas (t,R)RdR. We obtain the metal- 
licity at time t by its definition of Z(t) = M z (t)/M gas (t) 



Z(t)M„.Jt) 



(z M, 



out +f c 'R-*m z - e{ZM % - ezZIZ* ) dr 



(13) 
This is an integral-differential equation, which is easy to un- 
derstand in light of Z as an integral parameter with time. For 
a steady SF disk, its total mass is a constant, but metals are 
increasing. After some algebraic manipulations, we have its 
solution 



Z(f ) = Z max + (Zo - Z max ) exp 



rz 



(14) 



where Z max = (Z M out + f c K*mS)/(eiM m + e 2 Tl*), t z = 
M g3!i /{£\M m +E2R-*), and Z(0) = Z is assumed initially. For a 
time long enough, the last term will tend to zero, we have 



Z max =Z (1-A)+/ C (m z -Z () )A ; 



c f c m z 



Si + e 2 coM^ A 



-Ml 



(15) 



where A = TZ^/ieiM. + sjIZ*), M out = (l-f c )ll* +M. is used, 
M, is in units of M yr _1 , 7t* = c a M x m A , and c = 43.03 in- 
sensitive to M. for e = 0.5. The approximation is valid only 
for Zo -C m z . We assume for a comparison with observa- 
tional data that most of quasars reach the state with maximum 
metallicity, but the episodic age of quasars can be estimated 
from Equation (14) if the metallicity is measured accurately 
enough. 

Figure 4 shows a comparison of the model with data. Not 
only the metal-rich phenomenon in quasars is explained, but 
also the Z-L relation is reproduced by the present model. The 
large scatters of the Z-L relation imply different £1, £2, and c$ 
individually. 

3.3. Compact nuclear star clusters 

The present model may predict the formation of nuclear 
star clusters commonly found in local galaxies. They are 1-2 
orders of magnitude brighter than globular clusters (Cote et 
al. 2006), have extended star formation histories (Rossa et al. 
2006), complex morphologies (Seth et al. 2006), and follow 
the relation similar to the Magorrian relation (Ferrarese et al. 
2006). Though the SF disk discussed here is only on a scale 
of 1 pc, SF is ongoing inside the torus at a few 10 pc scale 
(Collin & Zahn 1999), which is invisible due to dust extinc- 
tion. Since AGN types are irrelevant to orientations of their 
host galaxies (e.g., Munoz Marin et al. 2007), SMBHs are 
randomly fed by the dusty torus. This is evidence for ran- 
dom accretion onto the SMBHs, which is derived from the 
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^-equation (Wang et al. 2009a) and is confirmed by calcu- 
lations based on the semi-analytical theory of mergers (Li et 
al. 2010). One stellar belt with random direction composed 
of old stars is left after one episodic activity accordingly. The 
appearance of the nuclear star cluster could be a natural con- 
sequence of the multiple episodes of SMBH activities (Wang 
et al. 2008; Wang et al. 2009a). The clusters are thus tightly 
related to the central SMBHs in the present model. Figure 3 
suggests that a mass ratio of SMBH and a nuclear star cluster 
is roughly 0.01-0.1 depending on e and j3 (i.e., star formation 
history), consistent with the observations (Seth et al. 2008). 
Detailed comparison with observed properties of the clusters, 
such as radial-dependent metallicity, SF history, will provide 
further clue to understand physics in galactic centers. 

All the emission from the SF disk will be thermalized 
with a typical temperature of Tsf ~ (Lsn/27t/?5 F <7sb) ~ 
2000 L S ^ 43 K in the near-infrared band, where osb is the 

Stefan-Boltzman constant, Lsn.43 = £sN/10 43 erg s" 1 is the 
SNexp luminosity. It then contributes a fraction to the ob- 
served infrared emissions in quasars, which is worth investi- 
gating further. 



4. CONCLUSIONS AND DISCUSSIONS 

We show that an SF disk supplied by the dusty torus is able 
to support an SS disk around SMBHs in AGNs and quasars. 
Such a model naturally produces the observed relationship be- 
tween metallicity and luminosity. As a natural consequence 
of the multiple episodes of activity through random accretion 
onto the holes, a nuclear star cluster will be formed from the 
remnants of stellar disks. The present model would provide a 
unified explanation of feeding SMBHs, production of metal- 
licity, and formation of a nuclear star cluster. 

We stress that the present paper deals with the stationary 
structure of SF disks switching to SS disks. A time-dependent 
model will give the evolution of AGNs driven by the SF disks 
for the complicated AGN-starburst connection. 
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